Surface energy and magneto-capacitance of superconductors under electric field bias 
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A superconducting layer exposed to a perpendicular electric field and a parallel magnetic field is 
considered within the Ginzburg-Landau (GL) approach. The GL equation is solved near the surface 
and the surface energy is calculated. The nucleation critical field of superconducting state at the 
surface depends on the magnetic and electric fields. Special consideration is paid to the induced 
magnetic-field effect caused by diamagnetic surface currents. The latter effect is strongly dependent 
on the thickness of the sample. The effective inverse capacitance determines the effective penetration 
depth. It is found that the capacitance exhibits a jump at the surface critical field. An experiment 
is suggested for determining the change in the effective capacitance of the layer. 

PACS numbers: 74.25.Op, 74.25. Nf, 74.20.De, 85.25.-j 



I. INTRODUCTION 



The gate voltage can be used to change the carrier den- 
sity of superconducting surfaces and therefore the critical 
temperature in the same manner as in field effect semi- 
conductor devices. This effect has been investigated for 
more than 40 years^ and has continuously gained exper- 
imental interest see also the overview^ The crit- 
ical temperature of thin superconducting layers can be 
controlled in this way by an electric field applied perpen- 
dicular to the layer i 1 ' 5 ' 6 ' 8 ' 9 ' High-T c superconductors are 
characterized by a low density of carriers such that this 
field effect is expected to be higher; this has caused a 
wide experimental activity ^ ^io, 11,12,13 g ucn field effect 
devices may be made even from organic and macromolec- 
ular films f 1 ^ 

According to the Anderson theorem the bias volt- 
age can change the critical temperature only indi- 
rectly via the electric field dependence of the material 
parameters i 4 ' 15 i 16 ' 17 i 18 ' 19 i 20 The influence of the electric 
field on the pairing mechanism is therefore to be ex- 
pected in the density of states for very pure and thin 
low-dimensional structure a 21 i 22 ' 23 analogously to the for- 
mation of sidebands in the density of states due to high 
fields^ Capacitance measurements on surfaces of high- 
temperature superconductors have revealed a so far un- 
known mechanism for electric-field penetration^ 5 - 

On the other hand, magneto-capacitance techniques 
arc used frequently to measure the influence of the mag- 
netic field 26 , e.g. to test spin-dependent electrochemi- 
cal potentials. Starting from the early reports on an in- 
crease of the capacitance^ 7 , for superconducting tunneling 
junctions 2 ^, the residual surface resistance of supercon- 
ducting resonators is still under discussio n 29 i 30 ' 31 since it 
becomes important for the question how short the elec- 
tron bunches can be in free electron lasers before the gen- 
erated wake fields disable the superconducting cavities^ 
Different mechanisms for such electron losses have been 
discussed in 31 . In this respect it is important to know the 



<8> 
B. 



FIG. 1: The slab superconductor of thickness d placed in two 
parallel capacitor plates at distance L. The superconductor 
extends infinitely in z direction. The magnetic field is parallel 
to the superconductor surface. 



explicit dependence on the magnetic field and the voltage 
bias. 

In this paper we investigate the magneto-capacitance 
in dependence on the magnetic field and the external 
bias by the electric field. We focus on magnetic fields 
around the surface critical field B c % since we expect that 
the external bias, which affects only the surface, has 
a relatively large effect on the surface superconductiv- 
ity. Most experimental activities are concentrated on 
the change of the surface critical magnetic field with 



temperature 



32,33,34,35,36,37 



In our study we suggest to 



consider these measurements under the influence of ex- 
ternal bias. To this end we will employ the Ginzburg- 
Landau (GL) equation with the DeGennes surface con- 
dition where only the latter condition depends on the 
external bias in agreement with the Anderson theorem. 

The experimental setup is shown in figure [T] The first 
electrode of a capacitor is the superconducting slab of 
thickness d. It is sandwiched between two plates of an 
ideal metal at distances L, which form the second elec- 
trodes. In real devices it will be necessary to consider 
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that the applied electric field affects also the end corners. 
Here we neglect the influence of the end corners assuming 
an infinite plane capacitor. Our aim is to evaluate the 
effective capacitance, the surface critical field and the 
surface energy in dependence on the applied voltage and 
the magnetic field. 

The total energy of the capacitor with area S is given 
by ^eoE 2 LS and an additional contribution coming from 
surface charges Sa. Then the inverse capacitance C of 
the slab is given by 
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The external magnetic field associated with the z axis 
and parallel to the superconducting surface is screened 
inside the superconductor by the diamagnetic current j. 
Thus a magnetic field profile is established perpendicular 
to the surface induced by the diamagnetic current. This 
induced diamagnetic current contributes to the inverse 
capacitance such that we have besides a genuine surface 
contribution C sur f also an induced part Ci n( j, 
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It will turn out that the induced capacitance is linearly 
dependent on the sample width d. 

Since the electric field penetrates the superconductor 
only near the surface it is of special interest to under- 
stand the surface superconductivity in the presence of an 
external electric field. The change of the upper critical 
field and the surface nucleation field has been calculated 
already for strong coupling. 38 The surface paraconduc- 
tivity and the change of the critical parameters due to 
an external field has been investigated, too i 39 ' 40 A shift 
of the critical temperature has been obtained 41 - due to a 
modified GL boundary condition and a variational solu- 
tion of the effective Schrodinger equation. In other words, 
the critical field B c \ is changed due to the change of the 
GL energy by the electric field42 Here we will investigate 
the surface energy problem of domain walls similarly and 
will employ the modified boundary condition to solve the 
GL equation variationally. We will obtain that the bulk 
critical field B c2 remains unchanged due to the applied 
electric field while the surface critical field B C 3 changes 
with the electric field. 

The paper is organized as follows. First we repeat the 
solution of the magnetic-field dependent GL equation un- 
der external bias by electric fields and calculate the sur- 
face critical field B c j, and its dependence on the external 
bias. With the help of the GL wave function the surface 
energy is calculated in chapter III and the effective ca- 
pacitance in chapter IV. Special attention is paid to the 
induced magnetic field effect in chapters III and IV. The 
selfconsistent treatment of the induced magnetic field is 
presented in appendix [D] which completes the proposed 
picture. We present fitting formulas for the magnetic and 



electric field dependence aimed for experimental verifica- 
tions. In chapter V we summarize and discuss possible 
experimental realizations. 



II. GL WAVE FUNCTION WITH EXTERNAL 
MAGNETIC AND ELECTRIC FIELDS 

An effective description of superconducting properties 
near the critical temperature is provided by the GL equa- 
tion for the wave function \P, 
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{-iHV - eA) 2 * + a* + /5|*| 2 * = 0, 



(3) 



which describes the ratio of the superconducting density 
to the total density n by l^l 2 = n s /2n. Here the mass 
is twice the electron mass, m = 2m e , and the charge 
is e = 2e e , the one of the Cooper pairs. If needed, the 
effective potential can be extended to lower temperatures 
than the critical one. 43 ! 44 

The GL equation is supplemented by the DeGennes 
surface conditions 4 ^ 
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where the extrapolation length b is sensitive to the treat- 
ment of the surfaces. This inverse extrapolation length 
1/b depends on the density of states at the surface, there- 
fore it is a function of the applied electric field E. In 
linear approximation, see appendix it reads 
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being of the order of few MeVs for conventional supercon- 
ductors. The dependence of l/(/?f c on the GL parameter 
k and the density derivative of the critical temperature 
suggests that the field effect is much larger for high-T c 
superconductors. 

The reduction factor T] is the ratio of the gap extrap- 
olated to the surface and the value at the surfaced Its 
value is of the order of unity and it is not essential for 
our discussion. 



A. Nucleation of the surface superconductivity 

At the surface critical field the superconductivity nu- 
cleates at the surface. Near the surface the effective wave 
function ^ is small and we can work with the linearized 
GL equation, omitting in (|3|) the cubic term, 
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with the boundary condition (|4|). 

We consider the geometry of a planar superconductor 
at d > x > as in figure [T] and assume a homogeneous 
applied magnetic field B a = (0,0, B a ). Since the system 
in Fig.[T]has translation invariance along the y direction, 
we use the Landau gauge of the form 



such that the wave function reads 



(8) 



The nucleation is possible if the parameter — a of ([7]) 
becomes equal to an eigenvalue e of the kinetic energy 
given by k^t (— «?iV — eA) 2 ip = eip. Since a changes with 
the temperature, a = a'(T — T c ), the eigenvalue e of the 
kinetic energy determines the nucleation temperature T* 
as T* — T c = -s/a'. 46 To avoid dual notation for the 
same quantity, we will treat the equation ([7]) as an eigen- 
value problem for a. Since a is negative, the nucleation 
temperature T* is always below the critical temperature 
T c in the absence of the magnetic field. 

Assuming the translation invariance along the y and z 
axes we can write the wave function as 



#(x,y,z) = tp(x) e iky e tqz . 



(9) 




Using ([9]) in the GL equation ([7]) we get a one- 
dimensional equation 

k — xj +q J ip + cttp = 0. 

(10) 

Any non-zero value of q results in the kinetic energy 
q 2 h 2 /2m which lowers the value of a making the nucle- 
ation temperature lower. The nucleation happens on the 
first possible occasion, i.e., at the highest allowed tem- 
perature. We thus take q = 0. 

The value of k determines the minimum of the 
parabolic potential and the eigenvalue a depends on the 
relative position of this minimum with respect to the sur- 
face of the slab. We have to find the wave vector k from 
the requirement of the highest nucleation temperature. 



B. Thick slab limit 

First, we assume that the superconducting slab is so 
thick that its surface superconductivity forms on both 
surfaces independently. In this case we can view the sam- 
ple as infinite and take the convergent solution into the 
bulk. We treat only the surface at x = 0. The surface 
at x = d is analogous. We should note that the profile 
for a general thickness has been solved by a calculation 
based on the Eilenberger equation for finite temperatures 
but without bias voltage4£ Here we restrict to a simpler 
approach but have included the bias voltage. 

It is advantageous to express the x-coordinate with the 
help of the dimensionless coordinate r 
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ip{x) = CD 9 (| + r ) 



(12) 



with the momentum To = — 2kl and the magnetic length 
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The parabolic cylinder function D v {t) solves the differ- 
ential equation 0, i.e4£ 



d 2 D v (r) 
dr 2 
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and the GL coherence length £ 2 = —h 2 /2ma. 

The boundary condition (fTS")) leads to a function v{tq). 
The maximal nucleation temperature is given by the 
maximal a — max[a] which is characterized by the mini- 
mum v = min[z/(To)] due to (|TB)) . Besides the obvious nu- 
merical search we can give directly a nonlinear equation 
for this desired minimum v'(tq) = 0. For this purpose we 
differentiate (|T5)) with respect to To, using the relations 
for the parabolic cylinder functions, D' v — tD u /2 — D u+ i 
and £)„+ 1 = tD u — vD v -\, to arrive at 




(17) 



With the solution v\tq\ of p7|) the momentum and cur- 
rent is determined due to To = — 2kl. In figure [2] the 
solution of (I17p is plotted. There is an asymmetry to 
be noticed with respect to positive, b > 0, and negative, 
b < 0, external bias. This will lead to very asymmetric 
curves in the surface energy later. 



C. Surface critical field 

The lowest eigenvalue v = min[^(T )] of ([H]) corre- 
sponds to the highest attainable critical magnetic field 



B c3 = max(S a ) = 
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where B C 2 is the upper critical field. The modified bound- 
ary condition does not influence B C 2 but makes it possible 
that — l/2<z><0 and a higher critical magnetic field 
B C 3 > B C 2 appears such that the superconductivity near 
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FIG. 2: The minimal Eigenvalue of the GL equation given by 
(|17|) as a function of external bias. 



the surface is enhanced in dependence on the external 
electric field. 

In figure [3] we present the result for the surface critical 
field (fT5|) versus the external bias We see that the 
external bias can enhance or decrease the surface criti- 
cal value depending on the field direction. Without ex- 
ternal bias the known GL solution B c3 /B C 2 — 1.69461 
is reproduced^ The strong coupling limit is somewhat 
larger resulting in the value 1.8i2£ We see that the elec- 
tric field can generate easily a value larger than 1.695. 
The experimental values compared with the GL theory 
and the theory of Hu and Korenman 33 are discussed in^S, 
which shows that the GL values are too small. 




condition ip 1 = 1/6 automatically. The constants a 
and To = — 2kl have to be determined by the minimal 
eigenvalue of (|7|). The normalized minimal mean eigen- 
value < a > can be obtained by minimizing the func- 
tional 
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with respect to tq and a. The resulting Z m 2 n = 2eB c3 /h 
yields 
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In the case of vanishing electric fields corresponding to 
the boundary condition 1/6 ps we recover the known 
results 51 , a = l/2£ 2 , x = l/%/2vra and B c3 /B c2 = 
\J ix jiji — 2) rs 1.66. The comparison of the expansion 
(|2U|) with the solution l|18p can be seen in figure [3] 



D. Variational wave function 



Let us now return to the full solution of (|T0|) . The wave 
function (fl"2)) specified by the value v = v describes the 
situation for the maximal nucleation temperature corre- 
sponding to the maximal magnetic field B a « B c3 . 




FIG. 3: The surface critical field B C 3 versus the external bias 
Q from linearized GL equation. The solution (|18|l (solid line) 
is compared with the expansion up to first (dotted line) and 
second (dashed line) order in 1/b according to (|20[) . 

To provide analytical expressions let us remind a sim- 
pler variational solution of the problem which was con- 
tributed by Kittel for absent external bias i 45 i 51 . We 
can extend this approach to external voltage bias and 
choose for the wave function the ad hoc ansatz ip « 
C exp (—a x 2 + x/b) which obeys the DeGennes boundary 



FIG. 4: The GL wave function of the condensate versus the 
external bias @. The normal GL wave function without ex- 
ternal bias is marked as thick line. The transition curve of 
£/l at B a = -B C 3 according to (|21[) is plotted as dotted line. 

In figured] we see how the external bias, 1/6, changes 
the wave function in dependence on the distance of the 
surface. The magnetic fields enters merely as a scaling of 
the spatial coordinate. For positive electric fields the su- 
perconducting density is diminished on the surface while 
for oppositely directed electric fields the surface super- 
conductivity is enhanced. 
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Since the wave function is strictly valid only near the 
transition line B a rs B c3 we have plotted in figure 2] also 
the transition line 



B a =B c3 




for x — £ as orientation. 

We will use the wave function (jT3J) as an ansatz for the 
variational calculation of the surface energy in such a way 
that the amplitude C serves as a variational parameter. 
This is motivated by the fact that near the surface the 
shape of the wave function is only slightly changed com- 
pared to the surface values but the amplitude decreases 
exponentially away from the surface. 

III. SURFACE ENERGY 
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FIG. 5: The surface energy in terms of the wall parameter 
l]28p versus the magnetic field and the external bias 1/6. The 
solid line denotes the transition curve B a = -B C 3 according to 
(|21[) where the surface energy vanishes. 



With the help of the GL wave function we can now cal- 
culate the surface energy. This surface energy is the inte- 
gral over the energy difference between the actual Gibbs 
free energy and the Gibbs free energy deep in the su- 
perconductor. Since the latter equals the one deep in 
the normal region when the field energy is subtracted, 
G(x — > -co) = G n o — B 2 /2^,o — G s o, we can write the 
surface energy as 



<r = J dx 
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= Jdx [G{x)-G s0 ] 
a 



dx 
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(22) 



The last counterterm is necessary in order to provide 
a consistent variational problem with the modified GL 
boundary condition (g]), for details see appendix [Bj 

The surface energy appears only if terms ~ |\I/| 4 are 
taken into account. Since the shape of the GL wave func- 
tion changes much less in terms of the applied magnetic 
field than the amplitude, we can now use our solution of 
the linearized equation, W = CD/}, to calculate the sur- 
face energy. For this purpose we determine the constant 
N, 



N- 



(3 



(23) 



such that 11221) takes a minimum. 



A. Limit of thick samples 

Assuming the limit of thick sample, d 3> I, we neglect 
in the first step the space profile of the magnetic field 



near the surface. This profile or induced field effect due 
to diamagnetic currents will be discussed separately in 
the next paragraph. In the same spirit, the upper in- 
tegration limit is taken as infinite here at this moment. 
The nonlinear Gibbs free energy (|22|) reads then 
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Here we have introduced 
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B = J drD%{r) 



and used a 2 / (3 = B 2 /hq. The minimum of 
N = A/B, 



(25) 



(26) 



leads to a surface energy a in terms of the condensation 
energy expressed in the critical field 
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with the wall parameter 



6(b) = -I 



B 
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In figure the wall parameter is plotted versus the 
external bias and the magnetic length. We see that 
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with increasing magnetic fields, i.e. decreasing magnetic 
length, the negative surface energy increases dependent 
on the external bias. Therefore the surface energy can 
be changed by the applied magnetic field as well as the 
external bias. The line of minimal eigenvalues of the 
GL equation (j2Tj) called transition curve is shown as well 
where B a = B c ^ and the surface energy vanishs. 

It is instructive to derive the wall parameter for the 
case without external bias, 
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with 
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FIG. 6: The surface wall parameter (|29|l for superconductor- 
vacuum boundary and zero external bias (|3ip versus the 
GL parameter (solid line) compared to the superconductor- 
normal boundary expression (|32|l of the literature (dotted 
line) . 
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J drDUr) j J drDt 



t (t) . (30) the superconductor- vacuum result (|3Tj) suggests a smaller 

value 



The upper critical field is related to the GL parameter 
B c2 = V2kB c and l 2 /f = B c2 /2B a . With the help of 
(HHJ) it is also easy to check that the surface energy (|29p is 
exactly zero for B a — B C 3 in the case of absent external 
bias. 

Assuming magnetic fields B a « B c in order to adapt to 
the situation of the superconductor-normal surface wall 
parameter we obtain 



lim 8(00) 

B„—>B C 



-^(0-41727- 



(31) 




An approximate treatment of the wall parameter of 
the superconductor-normal region is presented in ap- 
pendix [C] Neglecting induced diamagnetic currents 
the wall parameter S for both type-I and type-II 
superconductor a 45 ' 51 i 52 is approximated by (|C8[) 



(32) 



in terms of the London penetration depth A of the mag- 
netic field and the GL parameter k = A/£ with the co- 
herence length £. 

The result (|3ip for the superconductor-vacuum transi- 
tion can be compared with this superconductor-normal 
boundary ({52")) . In figure [H we plot (pITj) and (|52"j) and 
one sees the different places where the surface energy is 
vanishing. This vanishing of surface energy is connected 
with the transition from type-I to type-II superconduc- 
tivity. In the latter case the surface energy is negative 
indicating an unstable surface forming a vortex structure. 
While the superconductor-normal result (f3"2")l leads to 
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(34) 
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8^2/9 = 1.257 



(33) 



Without external magnetic field v — the 
superconductor-vacuum result (f3"Tj) coincides with 
the transition point between type-I and type-II super- 
conductivity. In other words the superconductor- vacuum 
boundary leads to smaller values of the transition be- 
tween type-I and type II than the superconductor-normal 
boundary with respect to the stability of the surfaces. 
We find that the applied magnetic field decreases the 
transition GL parameter. The type-II superconductivity 
extends towards values below k = l/\/2 = 0.7071. 



B. Finite width of samples 

When the applied magnetic field B a exceeds the upper 
critical field B c2 no bulk superconductivity is possible 
anymore but in a small surface region the surface super- 
conductivity occurs up to the surface critical field B c %. 
This critical field is dependent on the thickness of the 
sample 4 - . We have discussed the thick sample limit with 
respect to London's penetration depth, d 3> A, so far. 
Now we are going to consider the finite sample limit. 

Investigating a layer of finite thickness, d < 00 the 
Gibbs free energy (f2"4"]l with (f2"5"l) possesses an upper inte- 
gration limit d/l + TQ instead of 00. The superconducting 
surfaces are then also characterized by the appearance 
of diamagnetic currents. These diamagnetic currents in- 
duce magnetic fields, sec appendix ID] These fields con- 
tribute to the surface energy. We solve again the vari- 
ation problem of the free Gibbs energy (EM)) in order to 
obtain the optimum C. The expression for the surface 
energy (|22|) remains the same but the magnetic profile is 
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now spatial dependent, 

B(x) = B a +fi M(x) 



(35) 



The external magnetic field cancels in the difference of 
the Gibbs energies and is present only in the vector po- 
tential which takes now the form (ID6D, 



(36) 



AJx) = B a x + dx'M(x'), 



instead of 



As a result 
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Here we have used (|D13I) in the last line. Besides (|2"5)l 
we obtain now additional integrals, B' given by (|E7[) and 
V by (|E17|) calculated in appendix El 

The Gibbs free energy ([3T| contains now terms ~ C 6 
and the minimization yields a quadratic equation for C 2 . 
Instead of (|2"6")l we have 
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The resulting surface energy takes the form of (|27[) with 
the wall parameter 



A 2 



B + B 



-,giy) 



and 



3x 2 



3x - 1 + (1 - 2x) 3/2 



1 



x 
3 



(40) 
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From we see that the variational solution is only 
meaningful if y < 1/2. This specifies the lower limit on 
the thickness d m i n of our sample we can consider within 
this approximation. From ([39]) and with the help of (|E7|) 
and (lETTf we obtain 
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FIG. 7: The minimal thickness of the sample considered 
within the present approach versus the applied magnetic field 
and the external bias 1/6 in units of the coherence length £. 
The transition line according to (|21|l is given as well. 



where T given by JE9} and J by (|E18|) . This lower limit 
is plotted in figure [7] in terms of the magnetic length. 
One sees that we have practically a visible lower limit 
only for strong negative values of the external bias. 



1. Limit towards large sample thickness 

It is now interesting to discuss the limit of thick super- 
conductor probes. As shown in appendix [El both expres- 
sions B' and C are linearly dependent on the thickness 
d. We obtain for the wall parameter (|40|) 
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where T is given by (|E9[) . In the case of absent external 
bias one has 
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with £, 2 /l 2 — 2B a /B C 2, B C 2 — \[2kB c and the mean dis- 
tance 



d/l+To 
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J dTTD 2 (r)/ J drDUr). 
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If B a rj B c we have 
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-744326 -4 (0.41727 - n) 2 . (46) 
a 

We see that the surface energy vanishes at the same GL 
parameter as without induced fields (|3ip . 



IV. EFFECTIVE CAPACITANCE 

Now we return to the experimental setup of figure [1] 
and determine the expected contributions to the ef- 
fective capacitance ([T]). Besides the ideal capacitance, 
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FIG. 8: The induced part of the surface energy (|43|l in terms 
of the wall parameter (|28|l versus the magnetic field and the 
external bias 1 jb. The solid line denotes the transition curve 

B a = -Bc3 • 




f/b 



FIG. 9: The effective change of inverse capacitance in terms 
of (I48[) versus the magnetic field and the external bias 1 jb in 
units of the coherence length £. The transition line is plotted 
as in figure [8] 



Co = £qS/L we obtain contributions due to the external 
bias. 



S 

C cx 



e 2 dE 2 



Ln I d 2 



4(l-t) 2 = 4(l-t) 5 



Cox 



(47) 



where the temperature dependence t = T/T c arises from 
B 2 /2n = (e c /n)(l - t 2 ) 2 and k 2 (T) = 2k 2 /(l + t 2 ) in 
((6]) and we have scaled with respect to the temperature- 
independent coherence length £ = £(T)Vl — t 2 and (/? e i 
of ©. We abbreviate 



Lo 



5 2 



e 6Vel 2 ^o «eKg 



mc 
e 



3/2 



(48) 



Using the BCS expression for the GL parameter (3 — 
24ft 2 /(7C[3]nml.76 2 £! cs ) we can rewrite (HE) also into 
another form 



Lo 



t 77rC[3] 
1 6 



a B n 1.76k rj- 



dlnT c 



In i 



Ubcs (4g) 



with Bohr radius as and the Sommerfeld constant ckq. 

The contribution due to the external bias consists now 
of two contributions, 



1 

C C x 



1 



1 



c 



surf 



C 



ind 



(50) 



according to the thick sample limit of surface energy, (|2"T|) 
and (|28p . and the induced contribution due to diamag- 
netic currents, (|4"0)) . In the following we will discuss them 
separately. 



eigenvalue of the GL equation where the surface energy 
is vanishing. 

In order to provide easy to use formulas we can fit in 
terms of (05]) to obtain 



S 



C s , 



with 



/ 1 (a;) = 1.72(x- 0.39)e- 2 - 20 ^- - 39 ) 2 



(51) 



(52) 





e c /n 
[ M eV] 


Kg 


[10 28 m- 3 ] 


DbxT c 
dlnn 


tf In 7 
d ln n 


1/Vel 

1/MV 


Lo/eo 
nm 2 /pF 


Nb 
YBCO 


4.585 
750 


0.78 

55 


2.2 
0.5 


0.7422. 
-4.82 53 


0.4222 
-4.13 53 


4.52 
-207.5 


0.248 
2547 



TABLE I: Used material parameters 

Using as an estimate the parameters of table HI one 
finds for pure Nb L /e {t = 0) = 0.248nm 2 /pF while 
for YBCO one gets L /e (t = 0) = 2.547nm 2 /fF. This 
means that the expected change in the inverse capac- 
itance would be some nm 2 /pF. This should be com- 
pared to the trivial part of the inverse capacitance 
L/en. Measuring the distance L between the capaci- 
tor and the superconductor in mm one obtains L/e = 
112.9(L/mm)nm 2 /pF such that a relative precision of 10 3 
should be required to resolve experimentally the mea- 
sured effect. 



A. Limit of thick samples 

With the help of the surface energy shown in figure 
we can calculate the surface part of the inverse capaci- 
tance (1501) which is seen in figure O The larger the mag- 
netic field the larger is the inverse capacitance. For orien- 
tation we have plotted the transition line of the minimal 



B. Finite width of samples 

Now we focus on the induced effect by diamagnetic 
currents. The contribution to the effective capacitance 
([50)1 calculated from (|43|) is plotted in figure fTUl Com- 
pared to figure [9] one sees a different shape. Dependent 
on the size of the sample the induced effect can be larger 
than the thick sample limit. 
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FIG. 10: The induced part of the effective inverse capacitance 
in units of the inverse sample thickness scaled with k 2 £. The 
transition line is plotted as in figure [S] 



field. An experimental setup of such capacitance mea- 
surements is suggested supplemented by fit formulas. 
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The induced effect (l4"4l can be fitted analogously to 
ED as 



S 



Cind 



2; T / ; \ 3.25 



d £0 \£ 



"ind 1 ^ 



(53) 



with 



h ind (x) = -0.028c-°- 46 (*- ai9 > 2 

x (0.19 + 0.43x - 1.5a; 2 - 0.l6x 3 + .x 4 )(54) 

valid for £/b < 1.5. 

It is interesting to observe that both results, the ex- 
ternal (f2"B"]) as well as the induced one (|4"5)) show a nearly 
quadratic dependence 



C ~ Z 4 ~ B- 



(55) 



on the external magnetic field. This qualitative de- 
pendence on the magnetic field has been observed in 
magneto-capacitance measurements, 2 - 



V. SUMMARY AND CONCLUSION 

In this paper we have considered a superconducting 
layer under the influence of external magnetic fields par- 
allel to the surface and an electric field perpendicular to 
the surface. We found a nonlinear dependence of the 
surface critical magnetic field and the surface energy on 
the magnetic and electric fields. An effective capacitance 
is introduced which allows to measure these field effects. 
The diamagnetic currents induce a magnetic field profile. 
The selfconsistent equation for the induced magnetiza- 
tion in terms of a superconducting density profile is de- 
rived. The induced magnetization represents an impor- 
tant contribution to the surface energy and the effective 
capacitance. We report an explicit dependence of the 
surface energy and the effective capacitance on the layer 
thickness. The inversely-linear dependence on the thick- 
ness as well as the nearly quadratic dependence on the 
magnetic field of the inverse capacitance is in qualitative 
agreement with measurements of magneto-capacitance. 
We predict a similar behavior to be valid also in capaci- 
tor measurements on superconductors. The effective ca- 
pacitance is found to show a jump at the surface critical 



APPENDIX A: DERIVATION OF MODIFIED GL 
BOUNDARY CONDITIONS 

Let us outline the appearance of the modified bound- 
ary conditions (U|) . Therefore we start with the standard 
GL equation ^ 



— (-£W-eA) 2 # + a* + /3|#| 2 # = (Al) 

supplemented by the standard GL boundary condition 
that no current flows through the surface, 



Im 



(ihV + eA)§ 



(A2) 



The electric field will change the material parameters, 
i.e. the GL parameters a and j3. Linearizing the GL 
equation (|A1[) with respect to this induced electric field 
effect results in ^ = 'j/ + 5\E f where 'J obeys the GL equa- 
tion (|A1[) but with different boundary condition. To see 
this we consider the induced part of the wave function 
which is written as 6^/ — where ^ is proportional 
to the square of the Thomas-Fermi screening lengtbMis 
very small^ 4 - such that its explicit influence can be disre- 
garded. However it effects the boundary condition which 
translates now from (IA2I) for into the form for \P 



Im 



5A)* 



1 + * 



(A3) 



x=0 



Since ^ *~ E the external bias is changing the boundary 
condition of the GL equation into the form 



Im 



(iW - eA)W 



x=0 



1 

T 



E 

fie 



(A4) 



in linear order of the external bias. Since the vector 
potential is real in our case we obtain (0]). The ex- 
plicit form (U]) has been derived 4 ^ for layered super- 
conductors following the DeGennes theory 4 ^ or for bulk 
superconductors. 2 ? Please note that even without exter- 
nal bias the term l/i>o exists due to impurities and is 
usually small 4^ 
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APPENDIX B: SURFACE COUNTERTERM IN 
GL EQUATION 

For simplicity we consider the one-dimensional prob- 
lem with a superconductor in x > direction. The GL 
free energy density reads 

f{x) = f + a(x)V(x) 2 + ^'(x) 2 + |*(x) 4 (Bl) 

where a(x) = a + jp(tq + x/l) 2 and 7 = h 2 /2m. Search- 
ing for a minimum of 



and the domain wall parameter (|2T[) for B a w B c is 



5 / f(x)dx = 



(B2) 



with the ansatz ty(x) — CD u (x) with respect to the am- 
plitude C leads to (l23l) with ([25]) without counterterms 
and (1261). This minimum condition can be written also 



00 

= J (7*' 2 + a(x)t> 2 + /34' 4 ) 


00 

= y [-7*^" + a(x)* 2 +/3* 4 ] -7*(0)*'(0). (B3) 



This minimum condition differs from the Lagrange equa- 
tion of motion as the functional minimization of (IBID 



leading to the GL equation 

- 7*" + a(x)^ - /3<f 3 = 



(B4) 



(B5) 



just by the surface term \E f/ (0) 1 L'(0). For standard GL 
boundary conditions \P'(0) = this difference does not 
matter. With the modified GL boundary condition 
however, we have to add a counterterm to compensate 
this artifact, i.e. 



/(!)-> /(a) +7* , (0)*(0)*(a;). 



(B6) 



<5(oo) = Jdx 



r l + 2i-r'- + r- ( ^{x + T Q l) 2 -2 



(C2) 



One can simplify this expression by noting an addi- 
tional conservation law valid for one-dimensional GL 



equations^*. Multiplying the GL equation (|C1|) with 
ip(x) and integrating over x one gets exactly 

00 

'1 



dx 



e{^/i+r ) 2 ^ 2 +^ 2 )-^+^ 



= 0. (C3) 



Subtracting this relation from (|C2[) one has the simpler 
form 



5(oo) — J dx 
b 



-?P(x) 4 



M(x) 



(C4) 



The domain wall is characterized by an exponential de- 
cay of the magnetic field inside the superconductor and a 
decay of the superconductor wave function outside. One 
can calculate the domain wall parameter for two extreme 
cases analytically. 

In strong type-I superconductivity, k <C 1, no magnetic 
fields are inside, M[x) = —B a . Approximating £ 2 /I 2 ~ 
B a w and with the boundary condition at the surface 
"0(0) = and deep in the superconductor ■0(oo) = 1, the 
linearized equation J7]) has a first integral 



1 



£V(aO a = §(l-#s)) 2 



(C5) 



which can be easily verified by multiplying (|C1|) with 
2tj)'{x) and integrating from to 00. Taking this first 
integral (|C5|) into account one finds for (|C2[) 



(Si (00) = / dx [2£V 2 + (1 - ^ 2 ) 21 



OO 1 

2 J dx(l - tp 2 ) 2 = 2V2C J #(1 - V> 2 ) 2 



APPENDIX C: DOMAIN WALL SURFACE 
ENERGY 

Here we outline the derivation of the domain wall pa- 
rameter (|28|) without external bias b — ► 00 and for the 
situation of normal-superconducting boundaries as it is 
usually found in literaturei 45 ' 51 i 52 In terms of ^ the non- 
linearized equation ^ reads 







(Cl) 



= -V2£ = 1.89£ 



(C6) 



where we have used (|C5j) in the second last line once 
more. We see that in type-I superconductors the surface 
energy is proportional to the coherence length. 

The other case of extreme type-II superconductor is 
characterized by ip ~ 1 and an exponentially damped 
magnetic field profile with the London penetration depth 
A. Therefore (IC2D becomes 



$11(00) 



dx 



(C7) 
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and we see that the surface energy is proportional to the 
London penetration depth and becomes negative indicat- 
ing instability. We can combine (|C6|) and (|C7j) 



(C8) 



W2 
3k 



as an interpolation formula for both type-I and type-II 
superconductor o 45 ' 51 i 52 . 



APPENDIX D: INDUCED MAGNETIC FIELD 
DUE TO DIAMAGNETIC CURRENTS 

1. Selfconsistent magnetization 

We are going now to investigate the general form of 
the magnetization provided a profile of the wave function 
is given. The total induction becomes B(x,y) = B a + 
fioM(x,y) due to the induced magnetization (ioM(x,y) 
which is determined by the supercurrent j. For the sake 
of completeness we discuss also the y dependence. 

The external magnetic field and the magnetization are 
directed in z-direction 



{0,0, B a +fi M(x,y)) =rotA 



(Dl) 



from which one sees that A z = and the remaining pos- 
sible dependencies of the vector potential are A x {x,y) 
and A y (x, y) and 



dA v dA x „ , 

The current is given by 

j = (Kk-eA)-\1>(x,y)\ 3 
m 

1 „ fdM 8M n 
= _ ro tB= — ,- — ,0 

Mo \ oy dx 

from which we get the two equations 



(D2) 



(D3) 



= - e -mx,y)fA, 
oy m 

dM 

dx 



hk\ 



\^{x,y)\ 2 [A y --\. (D4) 



Eliminating the vector potential in (|D2[) with the help of 
(|D4[) we arrive at the differential equation for the mag- 
netization 

e ^r M = l (w^ M ) + 1 (wi M ) (D5) 

as the most general equation determining the magnetiza- 
tion profile. It represents a homogeneous linear differen- 
tial equation of second order which can be solved for the 
given geometry and a wave function ^(x 7 y) numerically. 



Here we restrict ourselves to a slightly simpler geome- 
try where we consider superconductors of large size in y- 
direction, see figure[T] This means that the dependence is 
only A v {x) and consequently M(x) which means A x = 
according to (|D4[) and further 



(D6) 



Now we introduce the dimensionless coordinates (jTTJ) and 



A y {x) = B a x + no dx'M(x') 



|^)| 2 = jjND 2 (r) 



(D7) 



where the wave function ip is given by the parabolic cylin- 
der functions D 2 {t) = D?(t) of IJT5J) only in the lin- 
earized ansatz while here D 2 stands for the square of 
the appropriately scaled general wave function. In these 
coordinates eq. (|D4[) leads then to the integral equation 



fM,M(x) _ N 



B c 
NB 

V5< 



T 

-J- [ cIt't 1 D 2 (t') 
V2kJ 



B r 



(D8) 



Differentiating (|D8|) twice and replacing one integral in 
the resulting equation by the expression obtained by dif- 
ferentiating (|D8[) once we arrive at the differential equa- 
tion 

h(r)y"(r) - h'(T)y'{r) + ah 2 (r)y(r) = 
2/(7-0) = 1, 2/(7-0) = ~aT h(T ) (D9) 



for the magnetization 

HqM{t) 

y(v = — 5 — + 1 

where we introduced the abbreviations 
h(r) = D 2 (r) 

NB r 



(D10) 



V2nB a 



(Dll) 



Please note that the trivial solution y — would mean 
total diamagnetism and complete Meissner effect. We see 
from the initial conditions in (|D9j) that this is ruled out. 
Instead we have a complicated profile. We can alterna- 
tively transform (|D9j) by z(x) = —y'{x)/{y{x)h{x)) also 
into a Riccati equation 



z'(x) — h(x)z 2 (x) — a = 



ar . 



(D12) 
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2. Approximate form 



Since our wave function (|12j) serves as a variational 
ansatz we can restrict the expansion here to terms up to 
the order C 4 . Therefore we obtain from (ID8D 



~ B, 

N 2 B. 



as 



TO 

r 

-JdT'Dl{r')Jdr'W , '(T^T'')Dl(T , ') + 0{N) 3 .(Bl3) 



This expansion can be alternatively considered as an ex- 
pansion for large applied fields B a /B c and/or large GL 
parameter k. 



APPENDIX E: INTEGRALS 



Here we introduce (|D13[) into ([57)) to obtain to addi- 
tional integral expressions beyond (|24|) proportional to 
B'N 2 and V'N 3 . In order to maintain legibility we in- 
troduce the integrals 



H(t) 
Q{r) 



T 

J dr'DUr') 

TO 
T 

J dr'r'Dlir 1 ) 

TO 
T 

| dr'r' 2 D 2 (r') 



(El) 



and write Gibb's free energy (|37[) with the abbreviation 



(E2) 



r>2 I at 

G = -l-S-\2NA-N 2 B 



2/z 
B C N 2 



W 



T 

J dT'T(j') 



2t / dr' 



2B a n 2 

1 

+ S j dT ' l:F {T"){H{T>)-H{T")) 



N 
V2k 



Hr') 



dr 



TO 

N 2 ^ n3b c r/ \ 
2k 2 v ; 2V2K 3 B a y 1 



x / dT'Hr')(H{T)-H{r')) 



B 2 

-I— c - (2NA - N 2 {B + B') + V'N 3 ) 



(E3) 



In the following the procedure will consist in rewriting 
the integrals into expressions which converge for w — > oo 
separating the terms explicitly dependent of w. This is 
achieved by systematically transforming the multidimen- 
sional integrals into integrals containing in each integra- 
tion a weight of D 2 . 

The additional term in (IE3j) proportional to N 2 reads 



(E4) 



W T 

2k 2 B' = J dr T 2 (t)+2tDI(t) J dr'Tir' 



Rewriting the first integral 



1 2 

T 2 {t) = J dr'r'D 2 ^) 



T T 

= 2 J dr'r'Dl(r') J dr"r" £>?(r") (E5) 

T T 

and interchanging integrations twice according to 

b t b b 

dr [ dr' = [ dr' [ dr (E6) 



a a 



we arrive? at 



w r 

2k 2 B' = 2 JdTT D 2 (t ) J <2tV-D?(t') (t-t'+w-t) 



wT 2 + 1 



with 



w 

1 = -2 f dTrDl(r)g(r) 



(E7) 



(E8) 



13 



where we have used (|E5[) once more and abbreviate Applying (|E6|) the first two integrals can be calculated 



H = H{w) 
T = T{w) 
Q = 9(w). 



(E9) 



The form (|E7[) shows that B' is proportional to the 
thickness of the sample since all remaining integrals con- 
verge to a finite value with C(exp (— w)). 

Next we calculate the term V'(V2n) 3 B a /B c = V' in 
(lE3l) with 



J drH(T) [2tT{t) - Q{t)\ = w 2 TU - TQ - wQU 

TO 

W 

+2 J <1ttDI{t)Q{t) (E14) 



and 



V' = / drD 2 (r) 



dT'T(T') 



T T 

h2r J cLt't 1 J dT"T{ T ")[H{T')-H{T")\ 

W T 

-2 J drT(T) J dT'Tir') [H(t) - H(t')} . (E10) 



J dTH(T)F(r) = T{vM - T) - QU + T 2 

TO 

W 

+2 J drDUr) [Q{t) - tT{t)] . (E15) 



In the first term we apply (|E5j) and interchange twice 
according to (|E6|) and in the second term we apply (|E6j) 
once to arrive at 



W T 



V =2 / dr J dr' {T{t)T{t') (H(w) - H(t)) 

TO T 

+F{™)Ht')[H{t)-H{t>)]}. (Ell) 



Employing (|E6|1 it is easy to see that 

dTUir) = wH(w) - T{w) 

J drJF(T) = w!F{w) — Q(w). 

TO 

Applying (|E6} and (fE12|) to (IeTT|) yields 

W 

V' = U(wT - Of - 2wT ( dTH{T)T{r) 



Successive application of (|E6j) and (|E12[1 yields after some 
straightforward but tedious steps 



JdTH(T)T{T) [tT(t)-G(t)} = j [wT{wT-2g) + g 2 ] 

TO 

W 

-\ J drD 2 {T)[g{r)-TT{T)] 2 . (E16) 

TO 

Using (|E14j) - (|E16j) in (|E13|1 we arrive finally at 



V 



(V2n)3B a 



wTJ -K, 



(E17) 



where 



(E12) 



W 

+2F J drH(T) [2tF{t) - Q( T )\ 

TO 
W 

-2 J dr Ti{r )J~{t ) [tT{t) - Q(t)] . (E13) 



W 

3 = T 2 -2 fdrDl{r)g(r) 

TO 

W 

K = - JdTD 2 {T)[TT{r)-g{T)Y 



TO 

-2T[Tg+l] 



(E18) 



We see that (|E17|) as (|E7[) are linearly proportional to 
the thickness of the sample. 
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